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ABSTRACT 

The  nonlinear  space-time  neutron  flux  equation  with 
negative  prompt  feedback  and  one-group  delayed  neutron  is 
reduced  by  the  use  of  a  nonlinear  transformation  to  a  par- 
tial differential  equation,  in  which  the  nonlinear  term 
represents  a  small  perturbation.   The  general  procedure  of 
solution  for  the  resulting  weakly  nonlinear  initial-boundary- 
value  problem  is  then  established  by  means  of  the  method  of 
successive  approximation.   Convergence  of  the  analytical 
solution  is  also  discussed.   The  solutions  to  a  slab  reactor 
core  and  a  cylindrical  reactor  core  are  investigated  here. 
Asymptotic  stable  equilibrium  states  are  derived  from  each 
of  these  solutions.   The  present  results  are  consistent  with 

those  obtained  from  previous  stability  analysis  for  the 

2 
generalized  buckling  K  greater  or  less  than  u  ,. the  lowest 

eigenvalue  of  the  associated  linear  Helmholtz  equation. 
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NOMENCLATURE 

r  position  vector 

x  dimensionless  position  vector 

t  time 

v  neutron  speed 

D  diffusion  coefficient 

k  multiplication  factor 

kc  multiplication  factor  of  an  infinite  reactor 

k  multiplication  factor  at  steady  state 

K  generalized  buckling 

C  concentration  of  one-group  precursor 

c  specific  heat  of  fuel 

P  * 

T  temperature 

T  temperature  of  coolant 

T  temperature  at  steady  state 

e  average  energy  released  per  fission 

h  convection  heat  transfer  coefficient 

C_  concentration  of  one-group  precursor  at  steady  state 

D  domain  of  a  reactor  core 

A  surface  area  of  a  fuel  element 
volume  of  a  fuel  element 


V- 


a       feedback  coefficient 
o 


B  fraction  of  one-group  delayed  neutron 

£  _  macroscopic  absorption  cross  section 

a. 

E  f  macroscopic  fission  cross  section 


T  dimensionless  time 

<f>  neutron  flux  distribution 

<(>  neutron  flux  at  steady  state 

A  decay  constant  of  one-group  precursor 

p  fuel  density 

ty  neutron  flux  difference,  <j)   -  <j> 

£  precursor  concentration  difference,  C   -   C 

6  temperature  difference,  T  -  T 

u,v  Eigenvalues 


ACKNOWLEDGEMENTS 

The  author  is  extremely  grateful  to  Dr.  Dong  Huu 
Nguyen,  Professor  of  the  Naval  Postgraduate  School,  for 
his  assistance,  encouragement  and  the  unselfish  giving  of 
his  time. 

Professor  David  Salinas  is  also  due  a  special  note  of 
thanks  for  his  comments  in  reading  the  rough  copy  of  the 
thesis  as  departmental  reader. 

The  author  is  deeply  indebted  to  his  wife  and  son  for 
a  long  time  that  they  offered  encouragement  and  exhibited 
exemplary  patience. 


I.   INTRODUCTION 

The  stability  of  nonlinear  reactor  kinetics  has  been 
investigated  by  several  authors.   Kastenberg  and  Chambre  [1] 
have  established  stability  theorems  for  the  nonlinear  space- 
time  reactor  kinetics  equations  by  the  method  of  comparison 
functions.   Garabedian  and  Lynch  [2]  and  Scalettar  [3]  have 
treated  the  nonlinear  equations  by  modal  expansions  in 
space  and  then  solved  the  resulting  ordinary  differential 
equations  numerically.   Kaplan,  Margolis  and  Harris  [4]  have 
applied  a  pertubation  technique  as  well  as  a  modal  expansion 
in  space  to  obtain  approximate  solutions  for  nonseparable 
transients.   Other  numerical  techniques  were  summarized  by 
Kaplan,  Henry,  Margolis  and  Taylor  [5]. 

In  this  thesis,  the  reduction  of  the  nonlinear  initial- 
boundary-value  problem  to  a  weakly  nonlinear  one  is  estab- 
lished by  a  nonlinear  transformation.   The  transformed 
equation  is  now  solved  analytically  by  means  of  the  method 
of  successive  approximations.   As  an  application  of  this 
method,  the  asymptotic  stable  equilibrium  states  are 
discussed  for  an  infinite  slab  reactor  core  and  a  finite 
cylindrical  core. 


II.   DERIVATION  OF  THE  SYSTEM  EQUATIONS 

A.   GENERAL  EQUATIONS 

In  this  thesis,  the  reactor  is  assumed  to  be  an  one- 
velocity,  bare,  homogeneous  and  Newtonian  cooling  system 
with  one  group  of  delayed  neutrons.   The  equations  describ- 
ing the  neutron  and  energy  distributions  and  the  concen- 
tration for  one-group  delayed  neutron  precursor  are  given 
by  Meghreblian  and  Holmes  [6]: 

±|i(r/t)  _  DVZ<P(C-t)  f  £a[(l-f#-']0(C't)  +  XC(Ct)       (2.1) 

£C(C/t)  =  -\C(r,-t)  +  p^>C,t)  (2.2) 

pS,|I(r,t)  =  eL^(C,t)  _  fc  A.  [T(r,t)-Tc(c,t)]  (2.3) 

where  the  resonance-escape  and  fast  nonleakage  probabilities 
to  thermal  are  assumed  to  be  unity.  The  boundary  conditions 
of  equations  (2.1),  (2.2)  and  (2.3)  are  given  by: 

.4>{$t)     =     C(£s,t)     =  T(nsi-t)  =  o   -for  -t>o  (2.4) 

and   the    initial    conditions    are 

<f>(r,o)  =-f(r);    C(r,0)   =  g(D    •    t(rfo)    =    <fccc)  (2.5) 

where  r   indicates  all  points  on  the  boundary  of  the  reactor, 

The  nonlinear  effects  on  the  kinetics  of  a  nuclear 
reactor  are  usually  described  by  a  single  parameter  ct  ,  the 
feedback  coefficient  of  the  reactor.   The  effects  of 


temperature  on  the  multiplication  factor  are  more  pronounced 
than  the  effects  of  temperature  on  the  other  reactor 
parameters.   Small  changes  in  temperature  cause  small  changes 
in  all  the  reactor  parameters.   However,  the  multiplication 
factor  is  very  close  to  unity  so  that  the  quantity, 
(l-3)k-l  is  very  sensitive  to  even  small  changes  in  k.   Hence, 
it  is  logical  to  assume  that  all  the  nonlinear  coupling 
appears  through  this  multiplication  factor  and  the  feedback 
equation  is  given  by 

.4(T)  =  k   +*o[T(rt)-T0(r)]  (2-6) 

The  dependence  of  the  multiplication  factor  on  the 
temperature  and  the  appearance  of  neutron  flux  <j>  in  equations 
(2.2)  and  (2.3)  make  the  system  defined  by  equations  (2.1), 
(2.2)  and  (2.3)  intrinsically  nonlinear  and  coupled. 

Introducing  the  new  dimensionless  variables 


T  =  v£  t 


*     =  \  /-§*  £ 


(2.7) 


the  above  system  of  equations  becomes : 

Sfyi/O-  Vf#(£/*)  +  [o-^k-iU(*/r)  +  AC(*,t)  (2.8) 


u. 


ac(x,T;  _  _  Jh_a*,~^  4-  £ik<K*<T)  <2-9) 

3T  vZa,  v 


W*  =  SrJ^}  ~  ^  [Tt*.*)-w]  <2-10> 
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Defining   the    constants 


«\  -  '"V  ;    22  -   pcvf     >       23  -     oTrrrEr?      >       A,  = 


P^X, 


f^^Z*      c_  v£. 


(2.11) 


equations  (2.8),  (2.9)  and  (2.10)  can  be  written  as  follows: 
2&&,*)      =  V^C*,?)  +  r(i-p)*,_l]#(*r'0  +  vAcC(i.-c)         (2.12) 


=  -^CCi,^)  +  y,4>C^/^) 


(2.13) 


JJt-*-)     =  ^fiL*,-*)  -  r3  [t(x(t)  -"£<*) 


(2.14) 


„     ..         ^.         ^ .    3d)    .   3C        ,  3T 
From  these  equations  ,  putting  •~-!-  =  0  ,  ~ —  =  0  and  -r —  =  0 

yields  the  steady  state  equations: 


0        -     7*0,(i)  +  [(i-f)to-l]^fi)+vA6,CodS) 


(2.15) 


0      =  -\c0(20  +.  y,  </>(<; 


(2.16) 


0      =  ^0eoO  -  ^[t0co-Tccm] 


(2.17) 


which  are  assumed  to  be  true  for  x  <  0. 

Subtracting  equation  (2.15),  (2. IP)  and  (2.17)  from 
equations  (2.12),  (2.13)  and  (2.14)  respectively  yields: 


§£l*.?)  =  V2^'^  +  [0-p)£-i]vK*/n  +  v7icE,(i,T) 


(2.18) 


St**) 

(7  L 


-A^t*,-?}  4-  Y,  pLZ,t  ) 


(2.19) 


9  L 


=  yavp(&,Tr)  _  #,#(*. -0 


(2.20) 
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where       ('-(O&e,-  '  ~  (i-p)-k.-\ 

and  M'U.t)   =  <f>^,t)~  4>0(&) 

}{*,t)      =   C0l,t)_CoC.s;  (2.21) 

and  k(T)  is  defined  by  the  feedback  equation  (2.6) 

,&(T)   =   &,0  -t-*o -£-(<, t)  (2.22) 

The  boundary  conditions  are  given  by: 

V(VO=  K**'T)   =■#(*->,?)  =0  (2.23) 

where  x   designates  the  position  vector  on  the  boundary  of 
the  reactor  and  the  initial  conditions  are: 

$(£,o)=   o  (2.24) 

Equations  (2.19)  and  (2.20)  with  initial  conditions  (2.24) 
can  be  solved  to  obtain: 

^)-K,   e'    ^C^f)^  (2.25) 

#(£/0  =  *zf  e"3f,^"'l')t  (*/•*')<**'  (2.26) 

It  is  customary  to  rewrite  equation  (2.26)  in  another  form: 
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•&(*,i)     =  fjH-r-t')^C&tx')J.x'  (2.27) 

where        3t(T-*')  =  ^  e  (2.28) 

The  explicit  form  of  the  kernel  (2.28)  depends  on  the  feed- 
back model  under  consideration  [7],   Substituting  equations 
(2.22),  (2.25)  and  (2.27)  into  equation  (2.18)  leads  to: 


1 


£-U,t;  =  Vz^(v(T)  +\(\-?)&o-i]y(i,^)  ^^-fit^C*^)  f3CCc-*W&,'T')dz' 

L  L  Jo 

+  vXe*JV'^'W,7<)^'       (2.29) 


which  describes  the  increase  in  neutron  flux  above  its 
steady  value  for  one-group  delayed  neutron  and  Neurizonian 
cooling  model. 

The  boundary  and  initial  conditions  for  equation  (2.29) 
are  given  by  equations  (2.23)  and  (2.24) 

vp(^,T)=   o  (2.30) 

q,(x,o)  =  ?(£)  (2.31) 

B.   ASSUMPTIONS 

Before  the  analysis  of  the  reactor  may  be  undertaken, 
two  further  assumptions  must  be  made  concerning  the  feedback 
model. 

The  first  assumption  states  that  the  reactor  temperature 
rises  instantaneously  with  the  flux  (i.e.,  the  prompt  feed- 
back model).   The  kernel  K(t_-c')   in  equation  (2.28)  becomes 
[7] 
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O-t(T-f)  =  rx   S(T-t')  (2.32) 

where  6(t)  is  the  Dirac  delta  function,  hence,  the  temper- 
ature difference  in  equation  (2.27)  can  be  written  as 
follows : 

-'o  o 

or  #(i/t)   =  ¥z^(±,~)  (2.33) 

The  second  assumption  states  that  the  temperature 
coefficient  of  a  reactor  is  negative.   This  is  the  case  of 
negative  feedback: 

|f  <  O  (2.34) 

Therefore,  the  feedback  coefficient  a   in  equation  (2.22) 

o 

is  also  negative: 

*0  <  0  -       (2.35) 

Inserting  equation  (2.32)  into  equation  (2.29)  provides: 

+■  vA^/V'^'W^')^'  (2.36) 

Let  fri_/i,WL_  l     =   K  (2.37) 


(2.38) 


vAt"*\       =    [3  (2.39) 
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Equation  (2.36)  yields: 


(2.40) 


and  its  initial-boundary  conditions  are  given  by  equation 
(2.30)  and  (2.31).   Equation  (2.40)  describes  the  increase 
of  neutron  flux  above  its  steady  value  for  a  negative  prompt 
feedback  and  one-group  delayed  neutron  model,  the  solution 
of  which  will  be  established  in  the  next  section. 
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III.   REDUCTION  OF  A  NONLINFAR  SYSTEM  TO  A  LINFAR  SYSTFM 

A.   NONLINEAR  TRANSFORMATION 

The  Wilhelra's  transformation  [8]  is  considered: 

*,*  ~n     -   eKaCS.-Q (3.1) 


or       U<4,-0     -    ag  ~  t£~\ r~  (3'2) 


Since  the  neutron  flux  difference  ^(x,x)  is  non-negative,  it 
is  logical  to  assume  that: 


W'CiS'O  >  °   fori  i.n>  2)    a-r,d       "E  >  o 
hence , 

i  -  £i  (i-  e*T)^(*,0  ^  o  (3.3) 

The  nonlinear  transformation  (3.1)  is  generally  valid,  since 
it  provides  a  unique  interrelation  between  tJj(x,t)  ^_   0  and 
u(x,t)  >_  0,  which  is  free  from  singularities  for  any  point 
(x,t).   Differentiating  equation  (3.1)  yields: 


1*  -  |^ir  +  ku.  _du.*)  §1. (3.4) 


Vzvp  -    e.  V  i>      ,   2*0- e  )e  (v^)  (3.5) 
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Substituting  equations  (3.1),  (3.4)  and  (3.5)  into  equation 
(2.40)  provides: 


\3l  » 


KT 


[-S(,-e")uJ 


KT_,Z 

e.    V  u. 

+ 

z^fi-e     )€    (vu) 

i-|(.-e    )u| 

<["i_  £(i-  e*  %)u 

i 

KT 

(3.6) 


Kt 


Dividing  both  sides  by 
give : 


[1  -  £  (l-eKT)u]2 


and  simplifying 


.-&(.-e  )< 


where 


<0     =   K  +  \.  >  o 


(3.7) 


(3.8) 


K>  o   and  7lc  >  c 


Let 


and 


7, 


(3.9a) 


(3.9b) 


Equation  (3.7)  becomes: 


3y  (<,-»-) 


=   vJuu^)+  i)x  (i-t)  +  ?,,  fi.r; 


(3.10) 
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The  initial  and  boundary  conditions  are 


"  (£..,  O   =   0 


(3.11) 


u 


f»,  <0   _   ^(*,o)   _  F(0 


(3.12) 


B.   COMPARISON  OF  TERMS  OF  THF  TRANSFORMED  EQUATION 

Before  solving  equation  (3.10)  by  successive  approxima- 
tions method,  it  is  necessary  to  impose  some  conditions  on 
the  terms  n,  and  n0  such  that  they  are  much  smaller  than 

the  other  terms  in  equation  (3.10).   A  comparison  of  the 

i~>t<\    K-c 
terms  of  equation  (3.6)  indicates  that  the  terms/ =p* J 6   and 

V  uje  '  have  the  same  order  of  magnitude  as  the  term  Kue 

Therefore,  the  nonlinear  terms  n,  and  n_  are  negligibly 

small  in  equation  (3.10)  for  any  time  0  <  t  <  °°  if  and  only 


if  the  ratios 


nl 

and 

n2 

KueKrp 

KueKrp 

are  small  compared  to 


unity,  i.e.,  for  x  in  D  and  0  <  t  < 


T?. 


Kue 


K"t 


iMteT-  i)(^u)' 


KleKT 


[i-*o-Ou]u 


« 


(3.13) 


and 


Kuar 


<Jrr-"')   u£*/?W 


1-sS  £,-«." 


'U 


(3.14) 


Since  u  and  Vu  are  regular,  as  Z    tends  to  zero,  the  numer- 
ators of  expressions  (3.13)  and  (3.14)  approach  to  zero, 
hence  , 
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11, 


Kue1^ 


->   o 


17, 


ku  e. 


For  the  time  being,  it  is  assumed  that  the  conditions  (3.13) 
and  (3.14)  are  satisfied.   They  will  be  verified  after  the 
solution  to  equation  (3.10)  is  obtained  in  the  next  section. 


C.   LINEAR  EQUATION  TO  BE  SOLVED 
Let  fe    =   1,  +  TJZ 

Equation  (3.10)  becomes: 


(3.15) 


—  V*u  (x,t)  +  g,[^  U(i,v>3        (3.16) 


u,(£s,t)  =   0 


(3.17) 


u(.*/°)  =  F(i) 


(3.18) 


The  transformed  equation  (3.16)  now  can  be  solved  by  the 
method  of  successive  approximations.   In  this  approach,  the 
small  nonlinear  term "g[t ,u (x,t) ]  is  treated  as  a  perturba- 
tion.  Hence,  the  initial-boundary-value  problem  is  written 
as  follows : 


|i!>    =   V*ul+  40,^,-0] 


U, 


I  (*;■*) 


(3.19) 
(3.20) 


";(*,<>)  =   F(i) 


(3.21) 
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where  ^, 

and 


£.  =  n.  +  TLs  -  2* tg"T-  ■> (?"i-f  +  fl*r, .. ■&  ( i - gKT)u,  l7e'^kaMlJ2'  (3. 22) 

thus,  the  reduction  of  the  nonlinear  initial-boundary-value 
problem,  equation  (2.40)  to  a  linear  one,  equation  (3.19)  is 
established  by  the  nonlinear  transformation,  provided  that 
the  conditions  (3.13)  and  (3.14)  are  satisfied. 
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IV.   ANALYSIS  OF  A  SLAB  REACTOR  CORE 

A.   SOLUTION 

Consider  an  infinite  homogeneous  slab  [Fig.  1] .   The 
medium  is  of  infinite  extent  in  both  y  and  z  directions  and 
of  width  2a  in  the  x  direction  (a  is  the  dimensionless 
extrapolated  half-thickness  of  the  slab) ;  the  origin  of  the 
x  axis  is  placed  at  the  center  of  the  slab.   These  specifi- 
cations reduce  the  spatial  dependence  of, the  flux  to 
variations  in  x  alone. 


-fcia. 


->  * 


Figure  1.   A  Bare  Slab  Reactor. 


Therefore,  equations  (3.19),  (3.20)  and  (3.21)  become 


ST  »  •   •  +  &i 


^(+8..?)  =  0 


(4.1) 


U- 


i(*,o)     =  F(X) 
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where  G>   =  o 


and 


-t  -  &«(«*-')  S^'P  (4.2) 


U  =  I,  Z,3,...) 


As  seen  in  Section  III,  equation  (4.1)  will  he  solved  by  the 
method  of  successive  approximations  where  the  zero  approx- 
imation corresponds  to"£,   =  0. 


th 
The  zeroth  approximation  u  (x,-c)  is  defined  by  the 


1.   Zeroth  Approximation  (£,   =  0) 


following  equations 


?U0  "^JJo 


u 


o 


(*,o)  =   Fc") 


Equation  (4.3)  is  of  the  homogeneous,  parabolic  type  that 
can  be  solved  by  either  separation  of  variables  or  using 
the  Green's  function  [9].   Throughout  this  work,  the  homo- 
geneous equations  will  be  solved  by  separation  of  variables 

Let         U*(*.r)       =   I(x)T(')  (4-4) 

hence,  %Mp        -       XT'       and  ^  =  X"  T        (4.5) 
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Substituting  equations  (4.5)  into  equation  (4.3)  and  then 
dividing  by  XT  yield: 


a 


X'   _  _X   =  _  uZ  (4.6) 

T   ~   X        " 

where  u  is  the  eigenvalue  to  be  determined  from  the  homo- 
geneous boundary  conditions : 

S(±a)  =  O  (4.7) 

Equations  (4.6)  can  be  written  as  2  ordinary  differential 
equations : 

X"  4-  ftZX   =   o  (4.8) 

T'  +.  j^T    =   o  (4.9) 

The  solution  to  equation  (4.8)  is 

X  (.*)        =    C,  cosyx  +  qsinpA  (4.10) 

where  c,  and  c„  are  2  arbitrary  constants.   The  boundary 
conditions  (4.7)  give 

C,c*5^g  +'  c^jrn^ta,   =   o  (4.11) 

and       C,  co^a-  --  £2s'nl*a-      =  o  (4.12) 

Subtracting  equation  (4.12)  from  (4.11)  leads  to 

This  implies  c_  =  0   or   sinya  =  0. 

If  sinya  =  0,  then  ya  =  n^,  n  =  0,1,2,...  and  equation 

(4.11)  provides 

C,  c»s n  it  =  o 
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or  C,   =   o 

and  the  general  solution  (4.10)  is  of  the  form: 

XCv)  =  Czs;n«_S  (4.13) 

This  solution  is  unacceptable  since  u  (x,-c)    is  non-negative 
for  -  a  <_  x  <_ci,    hence , 

C4   =   °  (4.14) 

and  equation  (4.11)  gives 

C,  co^^ua  —  o 

Since  c.  ^  0 ,  then 

Cojiutg.  =  o 

or  the  eigenvalues  are  defined  by 

U-n    =  (£n  +  i)l£      ,  fL  =  o,\  ,Z,..  .  (4.15) 

When  equations  (4.14)  and  (4.15)  are  inserted  into  equation 
(4.10)  the  solution  of  equation  (4.8)  is  obtained  as  follov/s 


X(x)   =  ^co^an  +  OJLx  (4.16) 

The  solution  to  equation  (4.9)  is 


T™   =  C^      **  (4.17) 


where  c,  is  an  arbitrary  constant.   Combining  equations 
(4.16)  and  (4.17)  gives  the  general  solution  to  equation  (.4,3) 
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u-„  (.*,-=) 


(4.18) 


where     are  the  eigenvalues  given  by  expression  (4.15). 

The  coefficients  B   are  determined  from  the  initial  condi- 

n 

tion:    u  (x,o)  =  F (x) 


(4.19) 


The  above  equation  indicates  that  the  B  's  should  be  the 
coefficients  of  the  Fourier's  series  expansion  of  function 
F(x)  ,  i.e.  , 


■J -a. 


(4.20) 


Replacing  B   in  equation  (4.18)  by  its  value  in  equation 
(4.20)  yields 


"  Cx.t)  = 


I 


«   -y  "£ 

n 
Since  the  series  E   e     cosy  x  cosy  £  is  uniformly  con- 

n  n      n 

n=0 

vergent  for  -  a  <_  x  <_  a  and  0  <  t  <_  °°  [cf .  appendix  A]  the 
interchange  of  signs  E  and  /  is  possible,  then, 


U0(x.~ 


&■$ 


(4.21) 


Let 


s(x'tiO  =  i£  e  ^co^ct^t, 


(4.22) 
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Equation  (4.21)  can  be  written  in  the  form: 


u  r  x,  -c) 


f  e(x/^,T)fr(3)a5 


(4.23) 


where  G(x,£,t)  is  usually  called  the  Green's  function  [8] 
associated  with  equation  (4.3) 
2 .   First  approximation 

The  first  approximation  is  the  solution  to  the 
following  equation: 


2t 


i^'  +  £,  fx.ti 


^x> 


u,(+a,T) 


(4.24) 


U.CX.o)  -   Fc«) 


where -£  ,  is  obtained  from  equation  (4.2)  by  making  i  =  1 
and  substituting  equation  (4.23)  for  u  (x,x).   The  solution 
of  equation  (4.24)  is  the  superposition  of  two  solutions, 
u   and  v  to  the  following  equations: 


(4.25) 


-3v 
9-r 


^+£/*,-r)  (4.26) 


UaLta.rf=     o 


W±a,-c)  =  o 


Uo(*,0)   =   F(x) 


u. 


u„  +  v 


v(*,  o^  =   o 


(4-27) 


Equation  (4.25)  is  identical  to  equation  (4.3),  hence 
the  solution  u  (x,t)  is  given  by  expressions  (4.18)  or  (4.23) 


26 


The    solution    to    equation     (4.26)     is    now   assumed    to   be: 

vO,t)      =     V    AaC?)  cos^ny  (4.28) 

the  y  ' s  are  determined  by  equation  (4.15).  The  source  term 

-£,  (x,t)  in  equation  (4.26)  is  then  expanded  in  terms  of  the 
eigenf unctions ,  cosy  x: 

i.Cx.t)  =  J  l^W***  (4.29) 

where         "&'  (T)   -  J  T^.  ft .  "t\cosu_£  d-E.  (4.30) 


Substitute  equations  (4.28)  and  (4.29)  into  equation  (4.26) 
to  obtain: 

oo   n  ***  op 

a=o  ri=o  ri:o 

dA 

where  A  (x)  =  -==- 
n       dt_ 


or 


\   [An(t)  +  ^An^)-^(T^Jco5ynx  =   o  (4.31) 


Since  the  eigenf unctions ,  cos  y  x   are  linearly  independent, 
equality  (4.31)  implies  that: 

Antt)4  V-n   An<-t)  =  "^1  tT)  (4.32) 

a.  —   o,  i  ,  2. , .  -  • 

The  initial  condition  is  given  by: 
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V(x,o)    =    O    =   Y^     A„(°)CoiKn 


which  is  equivalent  to: 


A^  Co)  =  O   ,   n  s  o ,  i ,  2., 


(4.33) 


^ 


Multiplying  equation  (4.32)  by  e  '   and  integrating  from  0 
to  t  yields  the  solution  to  this  equation  subject  to  the 
initial  condition  (4.33): 


where  ^  (t)  is  defined  by  expression  (4.30) 


(4.34) 


cLt 


or 


",T,.-Tj[*a" 


A.ct 


-Pnt*-*'. 


£,C£i*)co*f4*£<i*d;t 


(4.35) 


rv  =  0,1,1,5,.... 

Inserting  equation  (4.35)  into  expression  (4.28)  leads  to 
the  solution  of  equation  (4.26) 

Z — i  L     J0  J_a  J 

Since  £.(£,t)  is  well-behaved  for  regular  u(£,t),  the  uniform 

_y2(T-t) 

convergence  of  the  series   S   e        cosy  x  cosy  £ , 
a  n  n      n^ 

n=0 

[cf .  appendix  A]  allows  the  interchange  of  the  symbols  of 
summation  and  integration  in  the  above  equation: 
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•x   /■* 


Vtx.Tj  =  f     f    rAC>t)i|e"^(T*"*^xco5Fln5  i$<U     (4.36) 
-'o  ~'-a  *- 

From  the  definition  of  the  Green's  function,  equation  (4.22), 
equation  (4.36)  can  be  written  as 


(4.37) 


Jo     ■'-at 


As  seen  from  equation  (4.27),  the  first  approximation  is 
given  by 


U,C*,  ^)  =  u0c*,t}  + 


•Jo    -'-■a 


r-t)4u,tjd^t 


(4.38) 


3 .   i   Approximation 

It  is  noted  that  equation  (4.1)  for  the  i    approx- 
imation is  similar  to  equation  (4.24)  for  the  first  approx- 
imation with  the  same  initial  boundary  conditions. 
Furthermore,  the  Green's  function  associated  with  a  differ- 
ential equation  only  depends  on  the  form  of  that  equation 
and  the  boundary  conditions  [10].   Therefore,  the  solutions 
to  all  possible  equations  (4.1)  with  different  "£.  (x,x)  can 
be  deduced  from  expression  (4.38)  by  replacing-^.  (x,x)  for 
-L-L^t) 


u-l  (.*,T>   =  ^0(yt.-c)    -+ 


U   =   0,1,2, 


■Jo     -'-a 


(4.39) 


where  -£ .  (£  ,t)  is  given  by  equations  (4.2).   Once  u.(x,t)  is 
evaluated  from  equation  (4.39),  the  i    approximation  of 
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the  neutron  flux  difference  will  be  calculated  from  the  non- 
linear transformation  (3.1): 


KT, 


y.   #„  -c-)     _       e   u^  tV[j  (4.40) 


« 


Thus,  the  analytical  solution  for  an  infinite  slab  reactor 
core  with  negative  prompt  feedback  and  one-group  delayed 
neutron  is  found. 

At  this  point,  it  is  necessary  to  verify  the  conditions 
(3.13)  and  (3.14)  and  the  convergence  of  the  sequence 
{u.(x,t)}.   All  these  questions  are  discussed  in  the  next 
subsections . 

B.   VERIFICATION  OF  ASSUMPTIONS 

As  shown  in  part  III-B,  the  conditions  (3.13)  and  (3.14) 
are  justified  at  x  =  0  and  in  its  neighborhood  provided  that 
the  function  u(x,x)  and  its  derivatives  ^r—  are  continuous. 

oX 

These  conditions  will  be  verified  for  any  x. 

For  the  sake  of  simplicity,  the  zeroth  and  the  first 
approximations  are  used  here  to  calculate  the  values  of 
^(Xjt)  and  u(x,x)  respectively. 

Expression  (3.22)  yields  the  value  of  -£.  (x,x): 


•*,.=   <  +  1] 


where 


.2. 

r)'     =   -   fr*^-  '  lISiLi,  (4.41) 
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and 

7]'     fl*f,   -W.-p^u.  ]Z  f  £"C°  T"^arx,V)^/         (4.42) 
The  fundamental  mode  of  u  (x,x)  is  given  by  equation  (4.18) 

u0(x,t)   _    e>0  e   ccjs  M„  x 


lh 


and   so 


Let  a,(*)        -      2L&    cocia.x  (4.43) 

The  function  a, (x)  is  positively  small  for  -a  <  x  <  a  and  of 

a  -21 

order  of  —  B   since  a  is  very  small  (or  order  about  10    ) 
K   o  J 


O  <  -a,  7C    sL   60  «  I 

K 

Expression  (4.41)  and  (4.42)  are  written  as: 


(4.44) 


,     -**      _^     fK.^n^         ^ff)^/  (4.  46) 

^      =  f  Bee    [,-*,e   .*,*  jcos^J     *  ,,  (4-46) 

If   the   term   a,e       '      in   equations     (4.45)    and    (4.46)    is   neglec- 
ted  for   simplicity,    it   follows    that: 
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2.       J-  ^(-K-ZV^^ 


(4.47) 


^K<,T 


ri1  I     <   Pae       ,+a.cKje  £ ±E 


(4.48) 


a.      Case   K  =   y' 


Inequality    (4.47)    becomes 


t-nrr 


2.        4 


-|U,  T 


or 


r, 


Kue' 


KT 


^  K 


a^Qt 


-for    ar>^       T    y,    O 


Assume  a  -  10  21,   B   -  1010,   K  =  0.09,  X   =  10~6,  then 

'    o  c 


% 


Kue^ 


< 


-9 


x  /O 


or 


ti: 


K.ue' 


^   [   -(or  einy   t  >  <5 


In  the  same  manner,  expression  (4.48)  is  equivalent  to 


*   -k0t     2  r 


-It' 


(l  +  A,) 


^    -K0Tf    »  t' 


$  (^e^V-O    <      f- 


At 


32 


Replacing  K  by  (K+   )  in  the  above  inequality  yields 


and 


Mil    <   4re 


11: 


leu  a 


or 


1 


i<ue. 


>^T 


/       P  -for    a»^y 


T   >,    o 


Let    8*    -    10    8,      A      -    10~6    and   K   =    0.09,    then, 


I  kue*T      ^     io-' 


10     X  0.09 


I 

T 


t .  e- 


1 


<u  e. 


K.T 


[    //     I      -for   aM    T  j     o  <  t 


<    ^  (^.^o^ 


Hence,  inequalities  (4.49)  and  (4.50)  lead  to  the  following 
conclusion : 


Kui  e 


K.T 


<^V  I     -for  «li%3 


O  C  T  £•  <* 


(4.5-0 


This  means  that  the  nonlinear  transformation  (3.1)  may  be 
applicable  in  this  case. 


33 


b.   Case  K  <  y 

Inequality  (4.47)  can  be  rewritten  as 


Kj     <   *£*£*«■'** 


and  so, 


ti: 


ku  e. 


K.T 


</     ^fi  &„  f^e     '      <     IfL  &<>  \x\  ior  any  t  >, 


Let   y2   =    0.09,    K   *    0.08,    a    =    10~21,    B      *    1010,    X  10    6 

then, 


-II         to 
2x  lO      y  IQ   x  Q-Q9 

(o.od)1 


5.  x  \e> 


or 


o 


(4.52) 


Similarly,  inequality  (4.48)  turns  out  to  be 


*e'  =  -M{<2   fe     -  '     (4.53) 


Since ,  K   -  y 
o    o 


K  +  X   -  y 
c    Ko 


0.08  +  10~6  -  0.09  =  -0.01 


it  can  be  assumed  that 

(Ko-M-o  )   <   ° 

then,  inequality  (4.53)  yields: 
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or 


li 


KU  c 


ct 


< 


-(2.K-   ^*)t 


e 


K(^-'^ 


If  2K  -  11   >  0,  i.e.  ,  K  >  — ^—  ,  it  is  obtained  from  the 
o  —  —  2 

above  inequality  that: 


or 


hcwe' 


for   all   T   >    0 


< 


-8 


I  O 


-4 


o.  08(0.01")  8 


x  10 


This  implies  that 


*li 


xue 


Kt 


«   » 


for  all  T.  >  0 


(4.54) 


Expressions  (4.52)  and  (4.54)  prove  that  the  value  of 

1  -?  (x  t)  I 

j  ^   '|  can  be  considered  negligibly  small  compared  to 


,.   KT 

Kue  vS        2 

unity  for  any  0  <_  T  <^  °°  if  — =—  <_  K  <  y  . 


c.   Case  K  >  y 
o 

It  follows  from  inequality  (4.47)  that 


111!  <   z^t^  e 


for   any  t    >    0 


hence , 


n*. 


>  <  ti  e  K-c 


/      2o<  s0  H-f  s~  Mo       /        £<*  j.  K<? 


for   any   T   >    0 


-21  10         2 

Let  a  =   10         ,    B     =   10      ,    y      *    0.09,    K 

o  o 


0.093 


Then, 
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-21    io 


k  ae 


-  & 


or    LLl 


% 


sue 


KT 


« 


for  any  ■"  >  0         (4.55) 


To  prove  the  ratio   ™   is  very  small  compared  to  one  for 

Kue 

any  w   >_   0 ,  it  is  divided  in  three  different  cases  for  the 
sake  of  simplicity: 

I  .-     If   a,  €  <<   '      <-  ■  e 


-e  «    _l 


«*fi. 


or 


■X  «  7^7 


Expression  (4.48)  yields  the  inequality: 


hi  I    <  fh^ 


J 


/"•-f^"^,' 


or 


KoT     CK.-1-.^X 


y        a*  c  r    CK«-r»;x     "I   /     A* 

lil    <  Mi  e    0       -  'J<  £: 


.t. 


|C0-  ^ 


and 


Tlk 


ka  e 


^     /     £ e.     /       — £ —  for   any   x    >    0 


Assume  6*  =  10  8,  K  =  0.093,  K   -  y2  =  0.003 

'  '00 
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Then, 


Kue 


&._  If 


KT 


< 


10 


8 


-S 


a,  e> 


0.  0515*0. oo*> 

Ik-  ytf 


»        ' 


I        for    any    t    >^  0  (^.srg) 


<-  •  e 


7 


r 


It  is  obtained  from  inequality  (4.48)  that: 


1^1     <    ^,e       A,e 


i*-;^^' 


ca 


.e 


f  >  »» 


^e_^*-)T- 


<    f  £,  ^  e  r         )     ©  ^ 


K  A, 


tic-*p 


»M«  )t 


then, 


T 


Ku  &'" 


-lu^ 


K'  A, 


-for    ^iny    T    ;>, 


or 


It—     /     io  «lo   x  io       ^    , o     «   I 


for   any   m  (4.57) 


3._ 


if      a.,  e 


(K-tJ^)T 


*    i 


a.  •  & 


K- 


*r^(J<_)      -    7U7 


'o<8 
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then, 


hi  i  <  ^^ 


e      Jit  <  4(3  s0e   e 

CK^-jU.) 


J 


<     4£ 


8. 


-^«i 


(K.-H?) 


It  is  deduced  from  the  above  inequality  that: 


__i 


KK« 


ktt 


,.*.         -Kf 


/  ±  g         £         /  4ff for    an,     T    ^    o 


or 


i<ueM 


-6  -4 

/        ^*  [Q  x    io       <<    | 

"  o.  095*0.003 


(4.58) 


Therefore,  inequalities  (4.55),  (4.56),  (4.57)  and  (4.58) 
shows  that  the  conditions  (3.13)  and  (3.14)  are  satisfied 
for  any  t  >^  0 .   In  other  words,  for  all  cases,  it  can  be 
written  as  follows: 


t 


i<ue 


et 


«   I 


for  any  t  >  0      (4.59) 


C.   CONVERGENCE  OF  THE  SOLUTION 

Consider  the  sequence  of  successive  approximations 
{u. (x,t) | u. (x,t)  is  defined  by  equation  (4.39)} 


(4.60) 


(4.39) 


«C 


W(iar«-   G  (*,  £  /  f_  t)  =  -i-V 


_^C-r-fc) 


hn^  C°-T-n  I. 


cos  u  occoi 
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The  Green's  function  G(x,£,t)  is  continuous  in  -a  <_  x ,  £  <_  a 
for  0  <  t  <_  °°  [cf .  appendix  A]  and  bounded.  Thus,  the  con- 
vergence of  the  sequence  {u.  (x,t)}  reduces  to  that  of 

{"£■  (x,t)  }  for  0  <  x  <  °°. 

.    ti 
Since  the  ratio  ^—  is  very  small  compared  to  unity 

Ku.  e 

1 

for  0  <_  t  <_  °°  ,  we  have : 


K-C 


For  0  <_  t  <_  t,  ,  the  function  u.  (x,x)  is  analytically  well- 
behaved,  hence  bounded  by  a  number  M  and  the  inequality 
above  can  be  written  as 


«T. 


-£;    «    KM  <2 


which  implies  that  the~£[/s  should  be  bounded  by  a  quantity 


e' 


£.   <  l'  -for-    o  £  t     ±     T, 


In  addition,  equation  (4.39)  gives  the  following  recurrence 
formula: 

UUl-uL      =     f   Js^J.T-e^^tJ-ttCS.t)]^^    (4'61) 

a    —  £* 

U  =   0,1,2,-.- 

where  the'k.'s  are  evaluated  from  expression  (4.2).   Since 

_  o  1  _  p 

a  is  of  order  of  magnitude  about  10     and  3*  about  10 
[cf.  appendix  B] ,  the  first  term  in  the  right  hand  side  of 
equation  (4.2)  can  be  neglected  compared  to  the  second  term, 
i.e.  , 
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I.  e  fi*r,_  «.u.fx)u.f[L  j^Zl^^LlA-'  (4.62) 


It  can  be  seen  from  equation  (4.62)  that: 

since  the  integrand  under  the  integral  is  non-negative 
The  zeroth  approximation  yields  "£.   =  0  and  u  (x,t) 
Next,  equation  (4.61)  gives: 


U,  -  u  o  = 


-A 


where  G  (  *,  £  ,  t_  t-J   ^  o 

and  i,(^t)  >  O 

Therefore,  we  obtain 

or  ^,   }.   U0  (4.63) 

The  values  of  "c.  corresponding  to  u   and  u..  are  given  by 
equation  (4.62)  : 

2  /.i  -x^-r-t'; 


^  ^pT.-f^^^Ou^'f/ Uo^^'JaV  (4.64) 

The  values  of  u   and  u,  do  not  much  change  the  integrals 
in  equations  (4.64)  and  (4.65)  but  do  change  appreciably 


40 


the  factors  before  the  integrals.   Furthermore,  u,  >_  u  , 
hence, 

"£*  >  t,  (4.66) 

Combining  equation  (4.61)  and  inequality  (4.66)  leads  to 
the  following  inequality: 

or  u2  •%    u.,  (4.67) 

Following  the  same  procedure,  it  can  be  shown  that: 

and   u0  ^  u,  <.  ut  c  .  .  .  .  <  u.  <  lxUi  <  -  .  .  .  (4.69) 

The  sequence  {~Z,.}    is  bounded  and  monotonic  increasing, 
so  it  converges  [11] .   Therefore  the  convergence  of  the 
sequence  {u.}  is  proved. 

D.   ASYMPTOTIC  EQUILIBRIUM  STATES  FOR  AN  INFINITE  SLAB 
REACTOR 

Stability  is  one  of  fundamental  problems  in  nuclear 
reactor  design.   Some  questions  are  raised  here.   Is  the 
system  stable  if  the  steady-state  flux  <J>  (x)  is  perturbed 
locally  or  if  the  generalized  buckling  K  is  changed  by 
changes  in  the  absorption  cross  section  due  to  control-rod 
movement?   So,  this  sub-section  is  devoted  to  the  investi- 
gation of  all  asymptotic  equilibrium  states  (if  any)  from 
the  solution  obtained  in  paragraph  A  of  this  section. 
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1 .   Zeroth  Approximation 

The  zeroth  approximation  of  the  increase  of  neutron 
flux  above  its  steady  value  is  given  by  equations  (4.18) 
and  (4.40) : 

^(x.t)  =      4*>e      co^"*  (4.70) 


where  the  B  's  are  defined  by  equation  (4.20). 
a.   Case  K  =  y^  =  (^_} 

2 

All  (y   -  K)  are  positive  for  n  >  0  since 
n 

Therefore,  all  exponential  terms  except  the  term  in  (y  -  K) 
will  become  negligibly  small  after  sufficiently  long  time 
and  equation  (4.70)  reduces  to 

4/  (I   T)  >    B,eo^<     as  T  -5>  ^       (4.71) 

which  describes  the  zeroth  approximation  for  an  asymptotic 
equilibrium  state. 

For  the  case  of  no  delayed  neutron  (3=0)  and  no  feed- 
back (a  =  0),  expression  (4.71)  becomes: 

if>0(x,-t)  ^       B0C<5.s7u.x        3ls    T  — 5>  =^  (4.72) 

This  steady  state  is  exactly  the  same  as  that  obtained  in 
Reactor  Analysis  [6] . 
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2         2 
b.   Case  y   <  K  <  y. 

2 
All  exponents  (y   -  K)  with  n  >  0  are  positive, 

hence  the  corresponding  terms  in  equation  (4.100)  decrease 

2 
with  time  t  except  the  term  in  (y   -  K)  which  increases  with 

2 
time  since  y   -  K  is  negative.   Thus,  as  x  goes  to  infinity, 

ty    (x,t)  approaches  the  following  limit: 


^0Cx,t) 


B6-6.  cosujx 


-(Ho^-K)"" 


Since       21   B0g.~'"  coj/u0x        »    I  as     t 

the  stable  equilibrium  is 

^(x.-O  >   —  (4.73) 

For  the  linear  case,  the  reactor  would  be  called  supercriti- 
cal since  the  flux  i|j  (x,x)  will  tend  to  infinity  as  seen  in 
equation  (4.73)  by  putting  a  =  0.   Again,  this  result  is 

consistent  with  previous  Reactor  Analysis  [6] . 

2 
c.   Case  K'<  U. 

All  exponential  terms  in  equation  (4.70)  decrease 
to  zero  as  time  t  increases  to  infinity.   Thus,  the  increase 
of  neutron  flux  above  its  steady  value  is  zero  at  the 
asymptotic  equilibrium  state  for  the  nonlinear  case. 

vf^O.-c)  >       O         as   t  _=>   oc  (4.74) 

This  situation  is  called  subcritical  in  a  linear  reactor. 


43 


2 .   First  Approximation 

The  first  approximation  of  the  increase  of  neutron 
flux  above  its  steady  value  is  obtained  by  inserting 
equation  (4.38)  into  equation  (4.40)  and  putting  i  =  1: 


T   Q. 


—   (4-75) 


where  u  (x,x)  ,  G(x,£,T-t)  and  -£,  (£,t)  are  defined  by  equa- 
tions (4.18),  (4.22)  and  equation  (4.2)  respectively.   On 
the  other  hand,  -£,,  (£,t)  could  be  written  in  the  form: 

■£,(1,-0    =  Tt!(*..-fc)  +  nlu.o  (4.76) 

where  n,  (C/t)  and  nl(£,t)  are  determined  from  equations 
(4.41)  and  (4.42). 

Only  the  fundamental  mode  is  concerned  with  in  the 

following  calculations.   Therefore,  the  values  of  u  (£,t), 

3uo  ° 

v?—  and  G(x,£,T-t)  are  given  by: 

"of?,t)       =       8*.<2-        "J^,,*  (4.77) 

J-9(5,t)       =.   -B0^6e     si^p0^  (4.78) 

G(.>,i,T.(r)        =        4_e  Coi^Xtocp^  (4.79) 


Let  1,0,-0       =      e  <=>(x,S.JT-*)'fc|Urlt-;4$<Ur 

~* o  J. a 


(4.80) 


a 


<5C*.*,t-*  )£,  (S.fc^dtr  (4.81) 
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hence,  the  relationship  between  I, (x,t)  and  I_(x,t)  is 
established , 

I,(x,t)   =  eKXIzcx,-t)  (4.82) 

and  equation  (4.75)  becomes: 


VI^T)    =     AB(«,TUl,(>,t| 


,-^C«-a,e'r)uoci«.T)-^XlCK/t)+«illC}c#T 


(4.83) 


Substituting   equation    (4.77)    for   u    (x,x)    into   equation    (4.83) 
yields : 


^.Cv.t)  ^ corfi.x+X, (4.84) 


Taking  the  limit  of  the  above  equation,  as  t  tends  to 
infinity,  gives: 


&cg    ccs^+X, j  (4.85) 


It  is  apparent  that  the  temporal  behavior  of  iK  (x,t)  depends 
on  the  generalized  buckling  K.   The  limits  of  I,  and  I„,  as 
t  approaches  infinity,  are  evaluated  as  follows: 


C     C  f  *  C  "r 

Urn.     I  G(K.$,T-fcjfcr(*,e)«A$«lfc  =  jJiv^   j    \  c5(x,|,T-^^(S,t)dfc^ 

(4.86) 

Thus ,  the  problem  is  reduced  to  finding  the  value  of  the 
integral , 
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|    Gfif^^.tJ^  (J(t)clt    -for    T     lar^gc 
which   is    the   sum  of    two   integrals: 

6  O 

in  view  of  equation  (4.76). 

Equations  (4.41),  (4.77),  (4.78)  and  equation  (4.79) 
give : 


Let     A,(y,S)  =  -  1^.  S0  p*  c«cl«ox^i^o^sin2ju05 


(4.87) 


and  ct((f)       =        *1   80cosf>Lt£  (4.88) 


The    above    equation   can  be  written   in   the    form: 


f&ri'Jlfc    *     A,£*.|te.  e  ~    g: _       it  (4.89) 


Similarly,  taking  account  of  equations  (4.42),  (4.77),  (4.78) 
and  equation  (4.79)  yields: 
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Changing  the  order  of  integration  in  the  above  equation  and 
letting 

ft*(*iS)  =  A   &0   co^ju^x  co^zp.05  (4.90) 

leads  to  the  following  equation: 

J*i;-lt=A4e  J.,.4|e-,rw^«)Je  [,-«,e.W'   ]-£t*«  (4.91) 

a.   Case  K  =  M- 

As  t  is  sufficiently  large,  the  integrand  of  the 
integral  in  the  right  hand  side  of  equation  (4.89)  is  equiv- 
alent to: 

P.  -    <L  I 

~. ^t -tu«.  Kib   —    ,  (4.92) 

Substituting  equation  (4.92)  into  equation  (4.89)  gives: 


o  J,  l-f-  a 


dt_ 


or 


J  G^b    ol      A3(x,5)  e     T  (4.93) 


where       A5  (y,  s )  _.   Ai  (x,|)  (4.94) 

!  +  *,(*) 

In  the  same  manner,  equation  (4.91)  can  be  written  as 
follows : 
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_MNfV^-^X'  ^ 


|gi'  <it-  -   Ax  Cx,  g  )<2 


a  ^  °     ' '     ~"c 


e,  '  +  *,«) 


2.  (  H. S>0  * 


when  x  is  sufficiently  large   or  after  integrating: 

r^V^  a,A^rT,e^  _±V^    (4.95) 
But  equation  (3.8)  is  recalled: 


Ka  -  K +  ^c 


or  Xo-Ho2"  =  k-Ho^-4-^ 


since  K  -  p  ^ 


hence,         K.-jw^-    A  (4.96) 

Inserting  equation  (4.96)  into  equation  (4.95)  and  putting 


J\a(k,\)    =   Aa.Q-nM  (4.97) 


leads  to  the  following  result: 


or,  as  x  is  sufficiently  large: 

p^'^  ~       A+e"K°\  (4.98) 

Now,    adding   equation    (4.9  3)    to   equation    (4.98)    yields: 
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.t        -1-      «t  _UJ- 


(4.99) 


Integrating  equation  (4.99)  with  respect  to  £  from  -a  to  a 
will  give  the  integral  I„(x,t)  as  defined  by  equation  (4.81) 


.a-  /•t  ,,* 


PC  -U  "t 

I^Cx.-O  =1    I  &*,<{«:  4$  ^  A^C^)^  °^  (4.100) 


where 

5  *-'  - 

"'-•a. 


A5c>tJ  =  j[A5Cx.5)-*-A4C.<,gri^       (4.101) 


The  integral  I, (x,t)  is  evaluated  by  using  relation  (4.82) 
and  equation  (4.100) 

2 
Since  K  =  y  ,  then 
o 

I,Cx^  -  A^cxj-u  (4.102) 

Replace  equations  (4.100)  and  (4.102)  for  I„(x,x)  and 

2 
I, (x,t)  in  equation  (4.85)  with  K  =  y   to  obtain: 

.  L^      ^Cxf,)   =   Urn.   | &.«'*>*+*;  Ml | 

Since  the  term  A- (x) t  is  much  larger  than  the  other  terms 
in  the  right  hand  side  of  the  above  equation  for  x  large 
and  Aj.  (x)  is  always  positive  [cf.  appendix  B]  for  -a  <  x  <  a, 
the  limit  of  vp _,  (x,x)  becomes: 
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1<>yv   f  Ck,t)   _   J^  (4.103) 


•  ^ 


or 

V 

i.e.  ,  iK  (x,«>)  =  —  is  the  first  approximation  for  a  stable 

equilibrium  state  when  the  generalized  buckling  K  is  equal 

2 

to  the  lowest  eigenvalue,  y  ,  of  equation  (4.3). 

2 
b.   Case  K  >  y 
o 

2  2 

If  (K  -  y  )  is  positive,  so  also  (K   -  y  )  is , 
o  o    o 

since  K  =  K  +  X  . 
o        c 

For  t  large,  the  integrand  under  integral  (4.89) 
becomes : 

e   .  -  e- '      ~    £       =  _i_  (4.104) 

Similarly,  for  the  integrands  in  equation  (4.91) 


-CPo-k6)t'  -^o-k.Yt'  CK.-Ic)t' 

_€ I   .g 


i-^rtJe^^ttj^-"^    ^i)^'-K)T'  '  *^> 


_L  e*'**'  (4.105) 


and 


fl,U)e  (4.106) 


Substituting  equation  (4.104)  back  into  equation  (4.89) 
gives : 
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=     AcKs)e 


-r*: 


where 


(4.107) 


(4.108) 


Inserting   equations     (4.105)    and    (4.106)    into   equation    (4.91) 
yields    the    following   equation: 


GT}'   «LK     =    Aze 


*,(.*>    -<t 


then,  after  integrating  and  some  manipulations,  it  is 
obtained  that: 


:1I. 


(*-tf-*J 


K-  ft{ 


(4.109) 


The  exponential  terms  in  the  bracket  will  approach  zero  as 
T  tends  to  infinity.   Therefore,  equation  (4.109)  turns  out 
to  be : 


[G^i  Jib   ~ 


(K-i^)T 


o-  ,VC 


K  " 


V 


or 


1 


CK-Z^)T 


Let 


A7(x.£)  = 


AtC.«,TQ4.,Cg) 


(4.110) 


It  follows  that: 


I 


Gfl'A 


1' 


~  A?(k,  5)e 


(le-tp.)'* 


(4.111) 
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T 
Again,  the  value  of  integral  /   G&,dt  is  obtained  by  summing 

o 
up  equation  (4.107)  and  equation  (4.111): 


t_ 


f*-o  t  <■<-   ^  -> 


f  St/t   a  A^x,?)€_  °x  -f-  A7(».5)e 


2.   e        f^e   ^  T  +  A7 1  (4.112) 


[ 


When  t  is  sufficiently  large,  the  first  term  in  the  bracket 

is  negligibly  small,  compared  to  the  second  term  since 

2 

(K  -  u  )  is  positive,  hence 

Gt,lt         c       A?€.  (4.113) 

0 

The  value  of  I„(x,t)  is  deduced  from  integrating  equation 
(4.113)  with  respect  to  £  from  -a  to  a: 


rl     C         ,      ,  a     f<-2M«>T       (4  114) 


i 

a.  "  o 


A§(xl   =  J    A7(x/fe>)^5         (4.115) 


Where 


A_(x,t)  is  defined  by  equations  (4.110),  (4.90)  and  equa- 
tion (4.88). 

The  value  of  I.  (x,t)  is  determined  from  relation  (4.82) 

l.Cx.-c)   -  A&uie    n  (4.116) 

Inserting  the  values  of  I_(x,t)  and  I1 (x,t)  from  equations 
(4.114)  and  (4.116)  into  equation  (4.85)  yields: 
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tinv    t|/(  Cx,  -c  ) 


«.  -»pG 


Li, 


0oS  co^ju.y  ■<- A&<>oe 


HolT 


IK-  l"o   )T 


f     i  -  ft.cxie  '  "  ^  a,(«^e"  '  "'-  =iAg<>ie 


or 


Corn,    ip, ()*»"«) 


li  n\    , 


e [^60€  CoSfUX    j-  A  8  J 


-iK-^f,  )T 


/    -zfic-u2-)-!-         -Uie-pOt  -«-L/jn      .        _c-t  f 


Since  (K  -  y  )  is  positive,  all  exponential  terms  in 

the  bracket  go  to  zero  as  t  approaches  infinity.   On  the 

other  hand,  A0 (x)  is  always  positive  [cf.  appendix  B]  for 
o 

2 
-a  <  x  <  a,  thus,  the  limit  of  iK  (x,t)  for  K  >  y   is: 


c< 


(4.117) 


which  describes  the  asymptotic  equilibrium  state  of  the 

increase  of  neutron  flux,  c(>(x,t)  ~  4>  (x)  .   This"  is  exactly 

the  same  result  as  obtained  from  equation  (4.73)  for  the 

zeroth  approximation. 

2 
c.   Case  K  <  y 
o 

The  same  procedure  as  carried  out  in  parts  (a) 

and  (b) ,  the  stable  equilibrium  state,  will  be  found  by 

using  equations  (4.85),  (4.89)  and  equation  (4.91). 


-(fJ„L-KYt    -plb 

c     -  e. 


i  -  e. 


Kt 


I-  A,C 


rp?e 


(4.118) 


■+■  ^,  e 


r(H^K>t 


e*HS-wt_  -Kt+ 


a-. 
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The  above  expression  is  obtained  by  dividing  the  numerator 

-(Po-K)t 

and  the  denominator  of  the  left  hand  side  by  e 

2 
Since  y   -  K  >  0,  equation  (4.118)  becomes  for  t  large: 


i 


(po-<)t 


(4.119) 


eCf?-K>tr 


hence,  integral  (4.119)  can  be  evaluated  as  follows 


I 


Grj'  *Ur  cs.    A,  e     U       «^ 


A,e 


M/"-»< 


[,-e      J 


or 


pG^JU:   ^ 


a.  i).,  i>e 


rH 


2 


(4.120) 


where 


A  (>(i)  ,  AiC^X) 


(4.121) 


Furthermore,  integral  (4.91)  is  equivalent  to 


/ 


G*\' <U      ~  A2£      e  e       «Lt*t 


Ate 


rf» 


i-  ^  >— 


Cf^a-KolT 


It  is  assumed  that: 


H»  -  K°  > 


then,  for  x  large 
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—       Alo  ex,  *)  e 


6 

Since  ^       >      K0 

thus , 

J  ©£,  ^  -      A,0  ex.  Te,)e" 


—  <c.T 


and 


where  A„  fx)     _     I    AletK,?jJL5 


and 


(4.122) 


where         Alfl<K,"5)  =   A^x.g)  (4.123) 

(Mo4--^)2" 

Adding  equation  (4.120)  to  equation  (4.122)  yields: 

,2_ 


(4.124) 


J  3   z-"1"  I,,  T 

jcs£,^^  c=.   Al|Cx)^  '  (4.125) 


(4.126) 


I((x,Tr)  _  e  £xc*,-c)    =:  A,(c*)&      -.   A„t*;e  (4.127) 

Finally,  inserting  the  values  of  I„  (x,x)  and  I, (x,x)  from 
equation  (4.125)  and  (4.127)  into  equation  (4.85)  yields  the 
first  approximation  of  the  asymptotic  equilibrium  state: 
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I  I.  (  „      -CMo-hOT  _*«"* 

Li  rr\    ij/    ( x. ,  f  )      _       U no-  \  go  & 


coipj-x  -+-  A- ,,   cx)£ 


K     «         ~     T     K 

Since  all  exponential  terms  decrease  as  t  increases,  then: 

L^   ip>,  £x.,-r)  =   o  (4.128) 

Thus,  the  same  asymptotic  state  as  for  the  zeroth  approxima- 
tion is  obtained.   The  same  conclusion  can  be  made  even  if 

2  2 

K   is  assumed  to  be  greater  than  u  ,  instead  of  K   <  y  . 
o  '  o  o    o 

E.   PHYSICAL  INTERPRETATION  OF  RESULTS 

It  can  be  seen  from  the  zeroth  and  first  approximation 
that  higher  approximations  are  not  required  for  asymptotic 
equilibrium  states  in  a  nonlinear  system  with  negative 
prompt  feedback  and  one-group  delayed  neutron.   The  stable 
equilibrium  states  in  the  first  approximation  are  verified 
by  the  results  of  Kastenberg  and  Chambre's  stability 
analysis  [1] .   Furthermore,  in  the  absence  of  feedback,  the 
results  obtained  here  reduce  to  those  of  the  linear  reactor 
analysis  as  mentioned  in  part  D-l. 

The  equilibrium  states  are  established  in  a  nonlinear 
nuclear  reactor  at  different  levels  which  depend  upon  the 
value  of  the  generalized  buckling  K.   The  results  are 
interpreted  as  follows. 

For  t  <  0,  the  reactor  was  in  steady-state  with  flux 
4>0(x). 
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At  x  =  0,  the  flux  is  raised  to  an  amount  of  F (x)  from 
an  external  source.  What  happens  to  the  flux  distribution 
after  some  time  (t  >  0)? 

TT   2 

1.  If  K  >  (75—)  ,  the  neutron  flux  will  reach  the  stable 

~   A3. 

equilibrium  state  of  value  $(x)  =  <J>  (x)  +  —  after  the 
transient  given  by  equations  (4.39)  and  (4.40)  has  subsided. 

7T   2 

2.  If  K  <  (~— )  /  the  neutron  flux  will  return  to  the 

z  a 

initial  steady  value  (f>  (x)  as  equilibrium  state. 

On  the  other  hand,  the  equilibrium  state  is  weakly 
dependent  of  delayed  neutrons  since  the  fraction  of  delayed 
neutrons  (3)  is  very  small  compared  to  unity  in  the  expres- 
sion of  parameter  a: 


c<      =   I  ( I  -  (5  )<*<>  *z   I 


If  delayed  neutrons  are  neglected  ($*  =  0) ,  equation 
(2.40)  becomes 


?f  «.t)   =  *L±  (x,ti  -(-ic^c^t)-  oi  hAx-tO  (4.129) 

dH  "»** 

TT   2 

For  K  >  (^— )  ,  theorme  I  of  Kastenberg  and  Chambre  [1]  gives 
z  a 

the  asymptotic  stable  equilibrium  state  of  the  system 
(4.129),  which  is  the  positive  (nontrivial)  solution  of 
equation: 

(-irrv  f  i^jJi(*,T)  +  K+fx.t )-  PCtV^X/t)]   _   0  (4.130) 

K 
It  is  apparent  that  i|j  (x)  =  —  satisfies  this  equation,  thus 

is  a  positive  solution  to  equation  (4.130).   Therefore,  the 
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result  obtained  in  part  1)  is  verified  by  this  theorem  for 

7T   2 

K  >  (■=— )  .   Similarly,  theorem  II  of  Kastenberg  and  Chambre 

IT   2 

shows  that  the  equilibrium  state  for  K  <  (75—)   is  \p  (x)  =  0; 

£  a 

i.e.,  the  neutron  flux  returns  to  the  initial  steady  state 

cj>  (x)  .   This  is  again  consistent  with  the  result  in  part  2)  . 

2 

Theorem  III  states  that  if  K  =  y  ,  the  asymptotical 

stable  equilibrium  of  the  system  above  is  also  the  zero 
state  (i/j  (x)  =  0).   This  completely  differs  from  the  present 
result  as  described  in  part  1) .   From  the  physical  point  of 
view,  if  the  reactor  eigenvalue  is  not  changed,  then  the 

perturbed  flux  should  eventually  return  to  the  initial 

2 

equilibrium  state  4>  (x)  =  0.   The  fact  that  for  K  =  y  ', 

the  first  approximation  contains  a  secular  term  at  large 
r  [cf.  equation  (4.102]  indicates  that  this  is  a  special 
case  for  which  the  Wilhelm's  transformation  may  not  be  valid. 

To  draw  a  definite  conclusion  at  this  point,  it  appears  that 

2 
the  first  approximation  of  the  solution  for  the. case  K  =  y 

requires  further  considerations. 
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V.   ANALYSIS  OF  A  FINITE  CYLINDRICAL  REACTOR  CORE 

A.   SOLUTION 

Another  geometric  configuration  which  is  of  practical 
importance  in  Reactor  Design  is  the  cylinder  of  finite 
length.   Consider  therefore  a  bare  cylinder  of  dimensionless 
extrapolated  radius  R  and  height  21,    as  shown  in  Fig.  5.1, 
with  the  z  direction  along  the  axis  of  the  cylinder  and  the 
origin  of  the  coordinate  system  located  at  the  center.   The 
flux  in  this  reactor  depends  upon  r,z  (cylindrical  coordi- 
nates) and  t.   The  reactor  equation  is  now  [cf.  equation 
(2.40)] 


l±(,r,i,T)    =  V*>  +.  Kip-^t^^J 


'"  -   CT-T') 


->e< 


for 


0  £  r  <£   K. 


qj(r,i,-r')  J.T' 


(5.1) 


where 


Figure  5.1.   A  Bare  Cylindrical  Reactor. 


(5.2) 
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The  initial  and  boundary  conditions  are: 

(5.3) 
4^(o,  e,T)   lj-  -finite. 

4>  Cr,  +  l,t)    =       o 

Using  the  same  transformation  (3.1),  equations  (5.1)  and 
(5.3)  are  transformed  to  linear  equations  as  follows 
[cf.  equation  (3.31)]: 

2ii'Cr,8.,T^  =  V  utjcr,  e,T)  +  \   Cr,  fc,-c  )  (5.4) 

<-tiCr,*,©)   =  FCr,*) 

(5.5) 
Uc  (o,  e,T  )  ls   -f  ,'r^ite. 

u ;  (r,  ±  * ,  t  )  =  o 
where  -£e  =   o  (5.6) 

and 

The  methods  used  in  the  analysis  of  the  infinite  slab  are 
easily  extended  to  the  analysis  of  the  finite-cylinder 
reactor.   The  procedure  is  straight  forward  and  involves  the 
use  of  the  technique  of  successive  approximations  and  separa- 
tion of  variables. 
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1.   Zeroth  Approximation  (£   =0) 

o 

The  approuriate  form  of  the  transformed  equation 

which  gives  the  zeroth  approximation  is: 

1^     -     ±^(r^-°)     +  I*"0  (5.8) 


Uacr,  *.,   o  )   =  P(r,t) 
uo(£,*,t)   =   o 


(5.9) 


u0  (r;±  I  ,t)  -   o 


A  solution  of  this  equation  is  assumed  to  be  in  the  form: 

U„(r,ft,Tj   =  X(r)  Z  ca)T(r)  (5.10) 

The  substitution  of  this  expression  into  equation  (5.8) 
yields : 

T'<tn     l    a  l^\  Z"tt)  «.*  (5.11) 


^ 


T  (t  )    i"XCf 

where  V  is  a  constant  to  be  determined  from  the  boundary 
condition  (5.9)  or 

Z(~n  =  zu)  =  o  (5-12) 

It  follows  from  equation  (5.11)  that: 

Z"   +  |^£   =  O  (5.13) 
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and 

:l'       ^   i  I   /r2'\  _  _  ^  (5>14) 

T   r    rl  dr  {        J 

"jf  is  another  eigenvalue  to  be  determined  from  the  boundary- 
conditions  . 

X(0   =   o  {515) 

X  (  0  )   ii  -finite. 

Equation  (5.14)  can  be  written  in  2  different  equations: 

^.(rX')  +  ^rl   =   o  (5'16) 

and 

T'+C^+^T  =   o  (5.17) 

Solving  equations  (5.12)  and  (5.13)  gives: 

pt-n,  =   (Srn-i)TC  .  7^  =  0,1,2,. --  (5.18) 

and  .      2(4)  _   c,  ce*^*  (5.19) 

Equation  (5.16)  is  a  Bessel's  equation,  the  solution  of 
which  is  of  the  form: 

Xtr)   =   Cj0l^r)  +pY,(rr)  (5.20) 

where  C,D  are  2  constants 


X(o)   =   CJo  (o)  +  J>Y0(0) 
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hence,  D  =  0  since  y  (0)  =  °°  and  X(0)  should  be  finite 
[cf.  equation  (5.15)].   Furthermore,  the  boundary  condition 
X(R)  =  0,  yields: 


J0(^TR)   =  O 


(5.21) 


If  ■*$'   ,  y,  ,  /^-_  ,  .  .  rf   ,  .  .  .  are  the  positive  roots  of  equation 
(5.21),  the  solutions  to  equation  (5.16)  are  given  by  equa- 
tion (5.20) 


XCn)    =    CJ0  (tm  r) 


(5.22) 


and  the  solutions  to  equation  (5.17)  are 


T(r) 


E  £ 


"(<+  Hn  )X 


(5.23) 


where  E  is  any  constant. 

Combining  the  product  of  solutions  (5.19),  (5.22)  and 
(5.23)  leads  to  the  solution  of  equation  (5.8) 


u 


,Cr,*,T)  _£  A^ 


-(*m+f*«  ;** 


J0O^r)cosH_* 


(5.24) 


n\,n  s  o 


The  constants  A    are  defined  by  the  initial  condition 

m,n  J 


CO   _  o* 


U„(r.*,o)  =  F(r,*;   =   £  T  £  Am#  n  J0  CC^) 


n  =  o   m  i  « 


<^^£>M-»* 


^ 


This  equation  implies  that 


(  S-.  2.T) 
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Expression  (5.25)  is  known  as  the  Fourier-Bessel  series 

representation  of  F  (r) .   Since  V'  are  the  positive  roots  of 

equation  (5.21),  the  coefficients  A    should  have  the 
^  m,n 

form  [12] 


"  m,   n   = 


& 


f^jtanfO 


iJ.t^iK(i)Ari 


or 


A 


^VfflR) 


(   f^(t^M>~H.$^   (5'26) 


Substitution  of  expression  (5.26)  into  equation  (5.24)  gives 

the  solution  u  (r,z,x): 
o 


Kfoe/r)    = 


m.ixr  «  J    (^^)  ' 


Let 


G(ri*i*l.?,T)  = 


(5.28) 


nil"  =  • 


G(r  ,z,n.  ,£,x)  is  called  the  Green's  function  of  equation  (5.8) 
Then,  equation  (5.27)  becomes 


(5.29) 


2  .   First  Approximation 

The  first  approximation  satisfies  the  following 
equation: 
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2T 

— 

r 

W.cr, 

2, 

<0 

(5.31) 


ifsO4  ¥?'  -  *.cr..„  ,  (5.30, 

u  ,  (  o,  a,  t  )  is  f  int-te- 

where  "'t,  (r,z,T)  is  determined  from  equation  (5.7)  by  making 
i  =  1. 

The  solution  of  equation  (5.20)  is  the  sum  of  the  solu- 
tion u   to  equation  (5.8)  and  the  solution  v  to  equations: 

2v  _  J-?L(r™)    4-  ^X    +  1,  (5  32) 

v(r,a,  o  )    -   o 

/       N  •   r-   -L  (5.33) 

V(0,H.,  T  )      Li  finite 

v(ryii(-r)    =    O 

U,      =      Uo  +  V  (5>34) 

The  general  solution  to  equatoin  (5.32)  is  of  the  form: 
VM.-O  =  £  ^.J^  Jo^r)  ce'sr-2  (5.35) 

rh,nti> 
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Expanding^,  in  terms  of  J  (T  r)cosy  z  yields: 


va    »<■> 


^(^^'^ZZ  Cm,/-rlJ6(^r)c^u  4  (5.36) 


where   £   C„  J  (-^r)  are  the  Fourier's  coefficients  of  the 
_   m,n  o  m 
m=0 


above  series,  i.e., 


-e 


This  equation  implies  that  C    (~)  are  the  Bessel's  series 

m  /  n 

coefficients : 


C»,„tr>=        x\_  J    ^(t^H,?^)^^^^  (5-37) 

*R  J4  0'rnR)  J  o     -_./. 

Taking  derivatives  of  equation  (5.35)  gives: 


(5.38) 


m,  Tn  *  o 


Insert  equations  (5.38)  and  (5.36)  into  equation  (5.32)  to 
obtain: 
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<x> 


m  ,n  -  o 


This  equation  is  true  for  any  0  <_  r  <  R  and  -I    <    z    <    I   if 
and  only  if: 

4^   -*■  K^fO^   =   ^..-)         (5.39) 


The  initial  conditions  are  given  by: 


»r\.  r\  =  • 

or  &m„J0>>      =      °  (5.40) 

The   solution   to   equation    (5.39)    can   be  written   as: 
<W^    =     ^         £^^Mn)T      ^  +  Pm;     [C„,^)e  -Ur 

O 

where  D     are  constants  determined  from  the  initial  condi- 
m,n 

tions  (5.40) 

B^Co)   =  J>h..-  =  O 
Therefore,  it  is  deduced  from  the  above  equations  that 


Replacing  equation  (5.37)  for  C    (t)  leads  to 

III  f  XI 
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■»    r&  A 


e>    it) 


-rt-Vflfr-*' 


uXtt» 


*poft*))  2,(^.0*    '~ 'o^A^Lkr     (5.41) 


o      "  o       -£ 


Then,    it   follows    from  equations     (5.35)    and    (5.41)    that 

J.J.J,  «  fr_.  j/(ttr) 


•  -'«  "->i 


m.'X's  ° 


or,  in  terms  of  the  Green's  function  (5.28),  the  above 
expression  becomes: 


V(r,i,-r) 


f  f  j  i, h, 5,0^*,^  ,-*■-*)  «*£*»|^ 


(5.42) 


It  is  noted  that  this  form  of  solution  is  analogous  to  that 
one  derived  in  equation  (4.37)  for  a  slab  reactor. 

Now,  combining  equations  (5.29)  and  (5.42)  yields  the 
first  approximation: 


(5.43) 


U.Cni.-r  )          <-'0(r,*,T)    +     f 

J  t    - 

<S(r,  e,  /)  J,t-  t-)  i,  (n ,  ■§ ,  e)  i,\l*  dfc- 

3 .   ith  Approximation 

Since  equation  (5.4)  is  similar  to  equation  (5.30) 
with  the  same  initial  boundary  conditions,  the  ith  approx- 
imation which  is  the  solution  to  equation  (5.4),  can  be 
obtained  from  expression  (5.43)  by  substituting  -£..  (n  ,  E,  ,t) 
for  ^1(n,Cft)  : 
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n^  rt.   p*4 

Ui(r;*,Tj 

u0(r,*,t  )   + 

G(rV*,r)  .^.T-Owfi^.ty^^JU: 

J0  Jo  J  --e 

(5.44) 


Where  --£,.  (n,Cft)  is  defined  by  equation  (5.7)  .   Equations 
(5.44)  and  (5.7)  will  give  all  possible  successive  approxima- 
tions of  the  solution  to  equation  (5.4)  provided  that  -£,  =  0. 

Finally,  the  ith  approximation  of  the  increase  of  neutron 
flux  above  its  steady  value  is  calculated  from  the  transfor- 
mation (3.1) 


ICT" 

Vitr,£,t')    =  €  "«(r,g,T) (5#45) 

'-  £  (l-e*T)^(r,a-,T) 

Thus,  the  analytical  solution  for  a  finite-cylinder  reactor 
core  with  negative  prompt  feedback  and  one-group  delayed 
neutron  is  established.   The  convergence  of  the  series 
{u. (r,z,x)}  can  be  shown  by  the  same  method  as  indicated 
in  part  IV-C  for  a  slab  reactor. 

B.   DISCUSSION  OF  RESULTS 

From  the  solution  (5.44)  and  the  transformation  (5.45), 

the  asymptotic  stable  equilibrium  states  can  be  found  at 

different  levels  by  the  same  approach  as  applied  in  a  slab 

reactor.   The  proof  is  not  repeated  here  but  only  results 

are  summarized  as  follows. 

a.   If  K  =  y2  +-tf"2 
o    o 

After  a  long  time,  the  increase  in  neutron  flux 
above  its  steady  value  will  be  constant  with  time  and  reach 
the  value: 
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%(nt.x)    =  _M^^-!^ (5.46) 


where  A   is  given  by  equation  (5.26) 


Afl  = 


r£rl 


CR*  J,a(^0 


T)X(^r|JF(n/5)Cx.5^h5^^ 


-,1 


b.   If  K  >  y2  +  o/2 
o    o 

The  zeroth  approximation  for  asymptotic  equili- 
brium state  is: 

^,(r»  *  ,  «*  }     =  ■£■  (5.47) 

c  if  K  <  y2  +r2 
o    o 

%(r,£,*)         =   o  (5.48) 

It  means  that  after  the  transient  the  neutron  flux  will 
return  to  the  initial  steady  value  as  the  stable  equilibrium 
state.   It  is  noted  that  in  the  absence  of  delayed  neutron 
(3  =  0)  and  feedback  (a  =  0),  equation  (5.46)  becomes: 

H'o'tr,*,-)   =  A0ZU^)co^*  (5-49) 

for  K  =  y2  +  V2 . 
o     o 

This  result  is  exactly  the  same  as  that  one  obtained  in 
previous  work  [6] . 


2 .   First  Approximation 
a.   If  K  >  y2  +  T"2 

The  neutron  flux  will  approach  the  stable  equili- 


brium value: 
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9*(r,  *,oo)   _  6(r,i)    +  J£_  (5.50) 


c< 


as  time  goes  to  infinity.  - 

b.   If  K  <  y2  +  IT2 
o    o 


The  equilibrium  state  of  neutron  flux  is  the 
initial  steady  value  4>  (r,z) 

^(r,*,*,)   =   ^(r,0  (5.51) 

Therefore,  the  same  physical  interpretation  can  be  made 
for  a  finite  cylinder-reactor  core  as   for  an  infinite 
slab-reactor  core.   It  seems  that  the  asymptotic  states  for 
a  nonlinear  nuclear  reactor  do  not  depend  upon  geometric 
configurations.   This  is  true  since  the  transformation  is 
independent  of  coordinate  systems  and  the  transformed 
equation  is  maintained  in  the  same  form. 
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VI.   CONCLUSIONS  AND  RFCOMMENDATIONS 

The  objective  of  this  thesis  has  been  to  provide  an 
approach  for  directly  solving  the  nonlinear  space-time 
reactor  kinetics  equations.   The  theory  is  developed  for  an 
one-velocity,  bare,  homogeneous,  one-group  delayed  neutron 
and  Neutonian  cooling  system  where  the  temperature  rises 
instantaneously  with  the  neutron  flux  and  the  temperature 
coefficient  is  negative  (negative  prompt  feedback  model). 

The  fundamental  tool  in  this  analysis  is  the  Wilhelm's 
nonlinear  transformation  [8]  that  transforms  a  nonlinear 
equation  to  a  weakly  nonlinear  one  in  which  the  nonlinear 
term  represents  a  small  perturbation.   Then,  the  reduction 
of  the  later  equation  to  a  linear,  parabolic  type  is  estab- 
lished by  any  method  of  approximation.   The  successive 
approximations  technique  has  been  used  in  this  thesis  since 
it  yields  a  general  form  of  the  ith  approximation  for  the 
solution  and  it  is  easy  to  simulate  the  problem  by  digital 
computers  if  need  be. 

Furthermore,  this  transformation  is  independent  of  the 
system  of  coordinates.   In  other  words,  the  transformed 
equation  is  maintained  in  a  general  form  whatever  coordinates 
are  chosen.   Therefore,  the  same  conclusions  of  stability 
can  be  drawn  from  the  resulting  solutions  for  reactors  of 
different  geometric  configurations.   The  asymptotic  stable 
equilibrium  states  are  governed  by  the  relationship  of  the 
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2       2     2 
generalized  buckling  K  to  y   (or  y   +  -V  ) ,  the  lowest  eigen- 
value of  the  associated  linear  Helmholtz's  equation.   The 

V 

value,  <J>  (x)  +  — ,  is  the  stable  state  of  neutron  flux  when 
o      a 

2       2     2  2 

K  >  u   (or  y   +  Y   ) .   When  K  <  y  ,  the  neutron  flux  returns 
—  o      o    o  o 

to  its  initial  steady  value  tj>  (x)  after  some  transient. 
Then,  these  results  are  compared  with  those  derived  from 
Theorems  I,  II  and  III  of  Kastenberg  and  Chambre  [1]  in  the 
case  of  no-delayed  neutron.   As  mentioned  in  part  IV-E ,  the 

only  inconsistency  of  the  present  work  with  Kastenberg  and 

2 
Chambre' s  is  the  stable  equilibrium  state  for  K  =  y  .   But 

if  the  comparison  is  made  with  linear  reactor  analysis  (no 

delayed  neutron  and  no  feedback) ,  the  zeroth  approximation 

(cf.  see  part  IV-D-1)  will  give  an  asymptotic  stable  value 

which  matches  very  well  with  the  previous  result  [6].   So, 

for  future  work,  it  is  recommended  that  a  verification  of 

2 
the  first  approximation  at  K  =  y  '  by  another  approach,  say 

digital  computers,  be  made  to  insure  its  consistency  with 

the  theorem  III  of  Kastenberg  and  with  the  physical  point 

of  view. 

It  is  noted  that  the  successive  approximations  method 
can  be  applied  directly  to  equation  (2.40)  to  obtain  the 
solution  provided  that  the  coefficients  a  and  (3*  should  be 
very  small  compared  to  other  terms.   On  the  other  hand,  the 
Wilhelm's  transformation  reduces  equation  (2.40)  to  equation 
(3.10)  where  the  terms  n~  and  n_  become  smaller,  even  if  a 
and  3*  are  not  appreciably  small. 

Finally,  it  is  assumed  from  the  beginning  of  this  work, 
that  the  feedback  coefficient  a   is  negative.   May  the 
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transformation  work  with  a  positive  feedback  system  (a   >  0 

or  a  <  0)?   Moreover,  if  ty    >_   0,  u(x,t)  should  be  non-negative 

and  less  than  the  value,  -  =r= .   Otherwise,  the  func- 

a (e   -1) 

tion  ijj(x,t)  will  be  negative,  i.e.,  the  range  of  u(x,t)  is 
limited  in  this  case 


o    <  ul      <   _   __!< ,    c*  <   o 


and  the  singularity  may  occur  when  u(x,t)  is  equal  to 
K 


a(e   -1) 


at  any  point  (x,t).   Therefore,  it  is  not  con- 


venient to  apply  this  transformation  to  a  problem  with 
positive  feedback. 
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APPENDIX  A 

2 
UNIFORM  CONVERGENCE  OF  THE  SERIES   S   e     cosy  x  cosy  E 

n  n         n 

n=0 


Consider  the  above  series  where  y   =  (2n  +  1)  ■=—  .   We 

n  2a 

will   prove   that   this    series    converges    uniformly    for 
-a   £  x,    ?    <    a   and   0    <    t      ^_  t    f.  cc« 

y2TQ 
Expanding   e  into   Taylor's    series   yields: 


H^f. 


=         1+    rC^o   4-JL   ^T.)*^ 


hence,  e.  >     pft  ^ 


—  Fn  To  .  4-a. 

and  2.  ^        l 


Therefore,    inequality    (1)    can  be  written   as: 


C         CMf^xcec^l     <.         **  — '  -for    J 

ufc  i0        (.in  +  i;  [    T0  it  & 


It   is   known   that   the    series      E      ^  is    convergent,    then 

n=0  (2n+l)         2 

et     -y  x 

inequality  (2)  implies  that  the  series   S   e     cosy  x  cosy  E 

n=0 

should  converge  uniformly  for  -a  <_  x,  £  <_  a  and  x   <_  x  <_  °°. 
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Let       £(*.5,->)  =  ^  e.   uS^^»S[(^ 

n=  O 

Since  each  term  of  the  series  above  is  continuous,  the  func- 
tion S(x,£,t)  is  also  continuous  and  bounded  [11]  for 
-a  <_  x,  E,    <_   a  and  t   £  t  f.  oc- 
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APPENDIX  B 


1.   Calculation  of  the  Coefficient  A,,  (x)  in  Equation  (4.100) 


Assume  K  =  u.   =  {^—) 
i  o     za 


From  equation  (4.101); 


Ag-Cx)  =  J  A3  C«,  6)  <*£  -f-    A^  (x,  5)  4|, 


where 


K   «    '   ^ 

and 


a 
Let  us  evaluate  integral   /   A.,  (x,  £)  d£ 

-a 


r- 


-  .t!<.3tV.t»5tt.K[w^5l>,V'1  .ts   (3) 
Let      P   =   -Jl-  (4) 

Equation  (3)  becomes: 


77 


V*'2)^ 


4-So 


a. 


/V*f*-.x 


*V^ 


a., 


rp 


£ 


(5) 


The  integral  in  the  right  hand  side  of  equation  (5)  is 
separated  into  2  different  integrals  as  follows: 


I   -1 


=  iW^-pf'^f  V 


(6) 


But  the  first  integral  in  the  right  hand  side  of  equation 
(6)  can  be  calculated  from  C.R.C.  tables, 


J 


VT:rr  h  -  .% 


(7) 


The  second  integral  can  be  evaluated  after  changing  of 
variable: 


Let 


X.    =  y  +  p 


*lx  = 


^ 


hence , 


/»  P 


+< 


X 


9 


^Vlpy-tu  p2-     ^ 


(8) 


Since  it  is  obtained  from  any  integral  table  that 
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1      Jx.\l..^.  „*,, 


^  y^xS+.  7^>x.  +- 1  -  p1- 


(9) 


J     I/—  •X.V+.ip.x.-H- p*  ^  V^T        X 

sir.-'^.  (10) 


f  ^  -         1  cin-'f 


and       /  —  .     ..  =        '         cm-'f    P-^-^'-p 


t- 


=       '        .<■»-' />h-m  A  (ID 

Substitute   equation    (10)    and   equation    (11)    into   equation    (9) 
to   obtain : 

*V-x.\iyx.ir\-  p<-        j^    _^fT^  +  p5tw~V->/pvT*Siyr1^  "'-'  ")  (12) 

which   is    integral    (8) .      Thus   we   have 


(13) 


K 
Since  p  =   R   is  very  large,  hence, 


aB 
o 


■^  «  ' 


and  V^TT     ^O-^-V  a  p  -  _L  -  -J-s 


P       P      ^p3. 


Equation  (13)  can  be  written  as  follows 
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r'  

~  J^_ i_  (14) 

4p    5p*- 

Now,  inserting  equation  (7)  and  equation  (14)  into  equation 

(6)  yields: 

which  is  the  value  of  the  integral  in  the  right  hand  side  of 
equation  (5).   Thus, 


£A*0"'.0*3  =  -5r-^„«r. 


But  K0  =  2¥ 


Then 


r4 


J.>^^3   =  -i^t5-^-* 


(15) 


a 
Next,  evaluate  integral  /   A4(x,C)d^  where  A.(x,£)  is  given 

-a 
by  equation  (2) . 


=  £Te 


<J  CO  J 

<a.X, 
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Let  y  =  ^  =  iz  5 


dy  =  £  <*S 


hence , 


Cos   u  ^^4— 


it.  *■  *  ^ 


<=o 


(17) 


and  similarly, 


//-'V-JMj  =   *feJ^«W  =  .4iU) 


8   f. 
3   TO 


(18) 


Substituting  equations  (17)  and  (18)  into  equation  (16)  gives 


(19) 


Thus,  A5 (x)  is  obtained  by  summing  up  equations  (15)  and  (19) 


A^cx)  s 
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or 


V*-(-£+*^t  **-**.>.-■,... 


(20) 


Let  us  consider  the  order  of  magnitude  of  each  term  in  brac- 
kets.  To  do  this,  we  assume  the  following  data  [13] 


0*    =  5" 


5,  =  O ■  ooZ.  67crn' 

P     =     6.  06  «T7  §  / cm 


I8*   ° 


.  o«  3.5 


2*^   »•  * 


i?  e~2.  c»^i 


-/ 


ctj,  =   <3-  07?  dL/ QaK. 


Therefore,  we  obtain: 


K 


=  u„y-  -(%)*■«  °"i 


e<£ 


K  ci    -fe0_t      ^     ~6e    d     1.  0<jg£ 


T2=    _^ 


?f 


-it         2 

=       1,1x10  cru      4/c.     °K. 


V£*f<^ 


C<        - 


=    U.i4(jJ}*V>xl    ^  I  ^oY^  =     s~.r-^jp 


— -2.1 


cm     o«.c- 


7\      e      X 


0.  8  2-<+ 


_  t 


V<2T»  ^.  i^|>.\07  ^0  -oor' i- 


=       3.7.    *I0  C  <3-  i  rn  t*"  S 1  on  £/.*.!    ) 


-■Zo 


5     TCK 
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J$L    ^    £h     -    (!,-k0      =      o-ooic-x-  .  (21) 


-21- 


3  7TkXT   ^  -rr^J^Ac  J  /  (22) 


4-  X  J<-  -   /f.fex,o_7Jc^^o^^  (23) 

fi   =  (J-dsJ^Ac  =  o.oo  367.*  3 .2.  */ o  -  /.(x/o 

Thus,  comparing  equation  (21)  to  equations  (22)  and  (23) 
yields : 


-1  &il°      <za     it 

-3   7TK  Ac  Ac 

and        _L  JL  sL        <TC    £- 


and  the  value  of  A,,  (x)  is  approximately  equal  to 

/Vc(x)  .as.   ir  &-  c<^  3L  x  (24) 

Hence,  A,.(x)  is  finite  and  positive  for  -a  <  x  <  a. 

2.   Calculation  of  the  Coefficient  AR  (x)  in  Equation  (4.114) 

2 

Assume  K  >  y 
o 

From  equation  (4.115): 

4. 


A&^     =  J  A7(x,  50  4^ 


83 


where   A   =  .&z  ±i        =      Q>1 '<*   &°         coi«oX  co^C  H-o  5  ) 
7  *c(<-HS)         "■K^cLK.-p,?-)       ' 


Since  the  value  of  the  integral  in  the  right  hand  side  of 

expression  (25)  is  given  by  equation  (18),  hence, 

Aa(y)  _  A  (l*^g?-   oocX.                   (26) 

The  equation  above  shows  that  AR (x)  is  finite  and  positive 
for  -a  <  x  <  a. 
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